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ABSTRACT 

Context. The Crab nebula has been used as a celestial calibration source of the X-ray flux and spectral shape for many years by X-ray 
astronomy missions. However, the object is often too bright for current and future missions equipped with instruments with improved 
. | sensitivity. 

. Aims. We use G2 1.5-0.9 as a viable, fainter substitute to the Crab, which is another pulsar- wind nebula with a time-constant power- 

law spectrum with a flux of a few milli Crab in the X-ray band. Using this source, we conduct a cross-calibration study of the 
instruments onboard currently active observatories: Chandra ACIS, Suzaku XIS, Swift XRT, XMM-Newton EPIC (MOS and pn) for 
the soft-band, and INTEGRAL IBIS-ISGRI, RXTE PCA, and Suzaku HXD-PIN for the hard band. 

Methods. We extract spectra from all the instruments and fit them under the same astrophysical assumptions. We compare the spectral 
parameters of the G21.5-0.9 model: power-law photon index, H-equivalent column density of the interstellar photoelectric absorption, 
flux in the soft (2-8 keV) or hard (15-50 keV) energy band. 

Results. We identify the systematic differences in the best-fit parameter values unattributable to the statistical scatter of the data alone. 
We interpret these differences as due to residual cross-calibration problems. The differences can be as large as 20% and 9% for the 
soft-band flux and power-law index, respectively, and 46% for the hard-band flux. The results are plotted and tabulated as a useful 
' reference for future calibration and scientific studies using multiple missions. 

Key words, instrumentation: detectors — X-rays: individual (G21.5-0.9) 
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, 1. INTRODUCTION ences among various instruments. Benefits for the community 

• are two-fold: (1) The comparison result is useful in improving 

^ . From the successful launch to the end of the mission, all X-ray ±e calibration and in i nter p r eting data of other celestial objects 

7: ; observatories put tremendous effort into the flight calibrations obtained b diflferent instruments . (2) This stud will be a basis 

> . of their instruments using celestial objects to test their response tQ estsb]ish G2 1 .5-0.9 as a viable flight calibration source for the 

• -h ■ models and to monitor changes in their performance. Most such x . ray flux and the spectral shapei substituting the traditionally- 

^ ; efforts have been made by each observatory independently from used Crab nebula ( Kirsch et al 2005; Weisskopf et al 2010)j 

£ , the others, using similar methodologies and often the same tar- which h often t0Q brf ht for contemporary and future instru . 

gets. For the purpose of facilitating communication and compar- mentfJ wkh improved sen sitivity 
ison of results among various calibration teams, an international 

consortium for instrumental cross-calibration was established, Participating missions for the comparison studyare 

which is named the International Astronomical Consortium for the Chandra X " ra y Observatory dWeisskopf et al. | |2002j), 

High Energy Calibration (IACHECfl *f JJ™t^_G^^, Astrophysics Laboratory 

As part of the IACHEC activities, we present the results (INTEGRAL^^Winlder^et al. 2003) the Rossi X-«y Timing 

of a comparison study of G21.5-(X9 using the currently active 20 X P 7 ) rer s ( wi f^ eh r i a 3l t s e e t taL 2Q0 ^ ^ ^M^wSi 

X-ray missions. The purpose or this paper is to compare the 'i T 1 r i. n > — , ■ 

flux and the spectral shape in the X-ray bandpass of this time- Hansen et al|200lh observatories. 

constant power-law source and to identify the systematic differ- We use the following instrume nts: the Ad v anced CCD 

Imaging Spectrometer (ACIS; iGarmire et al.l l2003h on- 

* This work is based on the activity of the International Astronomical board Chandra, the INTE GRAL Soft Gamma-Ray Imager 

Consortium for High Energy Calibration (IACHEC). (ISGRI; [L ebrun et al. 2003) equipped with the IBIS telescope 

See http://www.iachec.org/index.html for detail. ( Ubertini et al] 120031) . the Proportional Counter Array (PCA; 
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iJahoda et all 1 19961) onboard the RXTE , the X-ray Imaging 
Spectrometer (XI S: iKovama et ap l2Q07h and the Hard X-ray 
Detector (HXD; iKokubun et all 120071: iTakahashi et all 120071) 
PIN component onb oard Suzaku, the X-Ray Telescope (XRT; 
iBurrows et alj 120041) onboard Swift, a nd the European P hoton 
Imaging Camera (EPIC) M OS-type dTurner et alj 1200 lh and 
pn-type (IStriider et all l2Q0lh CCDs onboard XMM-Newton. 
The dispersive X-ray spectrometers onboard Chandra and 
XMM-New ton are out of th e scope of this study. Readers can 
also refer to I S nowdenl (120021 ) . in which two past missions ASCA 
(GIS and SIS) and ROSAT (PSPC) are compared with Chandra 
(ACIS) and XMM-Newton (EPIC) using the same object. 

The instruments used in our study are divided into two 
groups: the soft-band instruments (Chandra ACIS, Suzaku XIS, 
Swift XRT, XMM-Newton EPIC-MOS and EPIC-pn) sensitive 
below ~10 keV and the hard-band instruments (INTEGRAL 
IBIS-ISGRI, RXTE PCA, and Suzaku HXD-PIN) sensitive 
above ~10 keV. The instruments with imaging capability are 
all the soft-band instruments with X-ray optics and INTEGRAL 
IBIS-ISGRI with coded mask technique. All the soft-band in- 
struments employ X-ray CCD devices, which have imaging and 
medium-resolution (E/AE ~ 30-50) spectroscopic capability. 

The plan of this paper is as follows. In §|2] we describe the 
X-ray properties of G2 1.5-0. 9 and discuss the advantages and 
limitations to use this source as a calibration source. In § [3] we 
present the basic properties of the participating instruments, the 
data sets, and the processing of the data. In §|4] we first set up the 
procedure for the event extraction and spectral fitting, which are 
common among all the instruments as much as possible (§ 14. U . 
We then present the fitting results of each individual instrument 
and justify the common procedure. We also discuss the effect 
of some major systematic uncertainties arising from the pro- 
cedure (§ I4.2l i. We finally assemble all spectra for joint fitting 
(§ !4.3l l. In §[3] we compare the results among all the participat- 
ing instruments and identify systematic cross-calibration uncer- 
tainties. We also compare the differences found in the present 
work with those in other studies. The paper concludes in §|6] 



2. OBJECT — G21 .5-0.9 — 

Pulsar wind nebulae (PWNe) are a small class of supernova 
remnants with intrinsically bright non-thermal X-ray emission, 
which includes the prototypical Crab nebula, 3C 58, Kes 75, PSR 
B0540-69, PSR B1509-58, and G21.5-0.9. PWNe are suited for 
X-ray flux calibrations for their constant total X-ray flux over 
human time scales. The brightness and the power-law index of 
the spectrum change spatially across a PWN, but the integrated 
spectrum can be represented by a single power-law in prac- 
tice. Therefore, they can also be used for calibrating the spec- 
tral shape, i.e., the power-law index in the X-ray photon spectra. 
Indeed, the Crab nebula has been used as a standard X-ray flux 
and spectral shape calibration source for decades. 

G2 1.5-0.9 is a PWN pow ered by the pulsa r PSRJ1833- 
1034 with a period of 61.8 ms (iGupta et al.ll2005b ICamilo et all 
120061) . The age and distance of the object ar e estimated as 
870 y r dBietenholz & Bartelll2008l) and 4.8 kpc dTian & Leahvl 
2008), respectively. The object is a bright X-ray emitter. The 
Einstein observatory resolv ed the object into several stru ctures 
of different spatial scales (iBecker & Szvmkowiak|[T98l1) . The 
EXOSAT and Ginga satellites showed t hat the X-ray spec- 
trum is represented by a power-law model (iDavelaar et all 1986: 
lAsaoka & Kovamalll990l) . arguing for the PWN nature of this 



G2 1.5-0. 9 has several advantages over other PWNe as a ce- 
lestial calibration source: (1) It is fainter than the Crab nebula 
by a factor of ~500, which places it in a reasonable flux range 
to match with the dynamic ranges of contemporary and future 
missions. (2) It is compact in size for its distance and youth, 
mitigating the effects caused by spatial differences of mirror and 
detector responses. (3) It has a very simple circular distribution 
of brightness and spectral hardness presumably because of the 
pole-on geometry, making the source and background extraction 
easy. (4) The spectrum is flatter than other PWNe, thus the ob- 
ject can be used to calibrate detectors over a wide range of en- 
ergies including those above 10 keV. (5) The soft (<1 keV) pho- 
tons contribute little to the spectrum due to a large interstellar 
extinction of >10 22 cirT 2 . Therefore, the calibration uncertainty 
stemming from the accumulation of contaminating material on 
CCDs, which plagues both ACIS and XIS, can be decoupled. (6) 
No pulsation is confirmed in the X-ray band, which indicates that 
the spectrum does not change in a pulse phase unlike the Crab 
nebula. (7) Because of these properties, G21 .5-0.9 has been used 
as a calibration source for Chandra, Swift, and XMM-Newton. 
We can exploit the wealth of existing data set. 

The source also has some limitations: (1) G21.5-0.9 is ex- 
tended by ~3' and is not adequate to calibrate dispersive spec- 
trometers. (2) The object is located close to the Galactic plane 
with a rich population of bright X-ray binaries, which makes 
source confusion non-negligible for instruments without imag- 
ing capability (PCA and HXD-PIN). (3) The spectrum suffers a 
large interstellar extinction, leaving few photons in the energy 
band below ~1 keV. Soft-band calibration sources need to be 
complemented, which is covered by othe r IACHEC studies (e.g., 
Beu ermannl20 06: Pluci nskv et al.ll2008l) . 
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Fig. 1. Chandra ACIS-S3 image of G2 1.5-0.9. The spatial com- 
ponents and the nearby source SS 397 are labeled. The source 
extraction region of 165" is indicated by the solid circle. 



Several spatial structures are known in the X-ray image. 
Figure Q] shows a Chandra ACIS-S3 view of the object. At the 
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center, the "core" is found (ISlane et alJl2000h . which is the neu- 
tron star with a power-law sp ectrum. The X-ray pulsation is yet 
to be found from t he pulsar dLa Pa lombara & Mereghetti 2002; 
ICamilo et al.l l2006>. The core is actually resolved with a size of 
~2", possibly rep resenting the pulsar wind termination shock 
(ISlane et al.ll2000l) . Surrounding this is the "plerion" with ~30" 
in size, which is a synchrotron nebula. The spectral hardness de- 
creases as the distance from the core increases, indicating the 
synchrotron loss of high energy electrons as a function of en- 
ergy. An even larger but much fainter structure is found around 
the plerion, which is called "halo" with -150" in size. The halo 
also shows a power -law spectrum with decreasing hardness as 
distance increases dWarwick et al.ll2Q0lT) . The origin of the halo 
is different from that of the plerion. One idea is that the halo is 
the du st scattered light from the plerion (Band iera & Bocchinol 
120041) . Upon these non-thermal features, a faint knot is found 
~1'.5 north of the core, which is called the "North Spur". Based 
on the thermal spectral fe atures, it is considered to be ejecta from 
the supernova explosion dBocchino et al.ll2Q05h . 

Some PWNe (e.g., the Crab nebula) show time-v ariable fil- 
aments. In G2 1.5-0. 9, dMatheson & Safi-Harbll20Toh have used 
archival Chandra data to detect some localized variable emis- 
sion near the pulsar. However, the contribution of these variable 
features to the total flux of G2 1.5-0. 9 is negligible of <1% and 
therefore our analysis is not affected by these features. 

An unrelated point source is found ~2'0 southwest of 
the cor e. The object is identified as an emission-line sta r 
SS 397 (Stephe nson & Sanduleakl \l9~TA IWarwick et all l200ll) . 
The source is a y Cas analogue, which is a small class of Be 
stars known to show weakly variable X-ray flux with extreme 
hardness. 

All these spatial structures can be separated only with tele- 
scopes with a good spatial resolution. Because the participat- 
ing instruments in our cross-calibration study include those with 
low spatial resolution or the lack thereof, we use a large enough 
source extraction region to encompass all of these structures. We 
ignore the contaminating thermal emission from the North Spur, 
which is very faint at a level of ~0.1% of the non-thermal emis- 
sion. We retrieved all the Chandra archive to examine the behav- 
ior of the unrelated neighboring source SS 397, and found that 
the source is variable by a factor of a few, but has a negligible 
flux even at the brightest level. The X-ray spectrum is dominated 
by the integration of spatially-varying power-law emission. We 
only treat the integrated spectrum, which we approximate as a 
single power-law spectrum. The reduced^- 2 values in our fitting 
will show that this approach is valid for the present data set. 



3. INSTRUMENTS, DATA, & PROCESSING 

3.1. Chandra ACIS 
3.1.1. Instruments 

The ACIS is an X-ray imaging-spectrometer consisting of the 
ACIS-I and ACIS-S CCD arrays. The imaging capability is un- 
precedented with a half-power diameter (HPD) of ~1" at the 
on-axis position. We use one of the CCD chips (ACIS-S3) in the 
ACIS-S array. The ACIS-S3 chip is back-side illuminated device 
sensitive at a 0.2-10 keV band. The chip has 1024x1024 pixels 
covering a 8'4x8'4 area. 




RA (J2000.0) 

Fig. 2. Layout of the field of views (solid squares) and the back- 
ground extraction regions (dashed rectangles) for each ACIS-S 3 
full-frame data set (CSO-6) in different colors. The source ex- 
traction region (solid circle) is displaced from the image center, 
leaving the other side of the image for the background extraction 
of a 5^5x1 '.9 region. 



3.1.2. Data 

The Chandra data were retrieved from the archive. A total of 
41 ACIS observations were available with G21.5-0.9 at various 
detector positions. We chose seven full-frame observations with 
the object close to but offset from the center of the ACIS-S3 chip 
(CS0-CS6; tableQ}, so that the chip encompasses the entire spa- 
tial structure of the source and leaves some area for background 
event extraction (figure|2|i. We also supplemented three ACIS-S3 
sub-array data (CS7-CS9; table[TJ to evaluate the effect of pile- 
up (§ 14.2. U . The data were taken in the timed exposure mode 
with a frame time of 3.2 s for the full frame (CS0-CS6) obser- 
vations or of 0.84 s for the sub-array (CS7-CS9) observations. 

3.1.3. Processing 

We processed the data based on the calibration database 
(CALDB) version 4.2 using the CIAO package version 4.2. 
The redistribution matrix functions (RMFs) and the auxiliary 
response files (ARFs) were generated for each data set. The 
RMFs take into account the spatial dependence of the detector 
response. The ARFs take into account the flux outside of the 
extraction region, although this is not a significant issue for the 
ACIS-S3 data as the resolved object is entirely contained in the 
source extraction region (figure [2]). 

3.2. INTEGRAL IBIS-ISGRI 
3.2.1. Instruments 

The IBIS telescope is a payload of INTEGRAL. It is composed 
of two layers of detectors. The upper layer is ISGRI sensitive in 
an energy range of 15-1000 keV band. ISGRI is equipped with 
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Table 1. List of data sets and their basic properties. 



Label 


Obser- 


Instru- 


ObsID 


Date 


'exp 


Band" 


c c 


Cnt" 








vatory 


ment 






(ks) 


(keV) 




bin 1 


(') 


(') 


CSO 


Chandra 


ACIS-S3 


1717 


2000-05-23 


7.5 


1.0-8.0 


28372 


50 






CS1 






2873 


2002-09-14 


9.8 


1.0-8.0 


36370 


50 






CS2 






3700 


2003-11-09 


9.5 


1.0-8.0 


34749 


50 






CS3 






5166 


2004-03-14 


10 


1.0-8.0 


35777 


50 






CS4 






5159 


2004-10-27 


9.8 


1.0-8.0 


36729 


50 






CS5 






6071 


2005-02-26 


9.8 


1.0-8.0 


34968 


50 






CS6 






6741 


2006-02-22 


9.8 


1.0-8.0 


35947 


50 






CS7 






1553 


2001-03-18 


9.7 


1.0-8.0 




50 






CS8 






1554 


r\ f\f\ 1 nn /i| 

2001-07-21 


9.1 


1.0-8.0 




50 






CS9 






3693 


2006-05-16 


9.8 


1.0-8.0 




50 






ISO 


INTEGRAL 


ISGRI 




2003-2008 


3132 


18-150 


3.39xl0 6 








RPO 


RXTE 


PCA 


20259-01-01-000 


1997-11-08 


21 


5.0-30 


1. 27x10" 








SIO 


Suzaku 


XISO 


104023010 


2009-10-10 


40 


1.0-8.0 


81297 


100 


5.0 


7.0 


SI1 




XIS1 




2009-10-10 


40 


1.0-8.0 


89436 


100 


5.0 


7.0 


SB 




XIS3 




2009-10-10 


40 


1.0-8.0 


83310 


100 


5.0 


7.0 


SPO 




PIN 




2009-10-10 


30 


15-70 


12684 


800 






sxo 


Swift 


XRT 


00053600001 


2006-08-13 


17 


1.0-8.0 


10593 


20 


5.0 


7.0 








00053600002 


2006-08-15 




















00053601001 


2006-08-23 




















00053601002 


2006-08-24 














SX1 






00053600004 


2007-05-09 


26 


1.0-8.0 


15767 


20 


5.0 


7.0 








00053600006 


2007-05-11 




















00053600007 


2007-05-16 




















00053600008 


2007-05-17 




















00053600009 


2007-05-29 




















00053600010 


2007-05-31 




















00053600011 


2007-07-04 




















00053600012 


2007-06-28 














SX2 






00053600021 


2007-10-06 


28 


1.0-8.0 


17999 


20 


5.0 


7.0 








00053600025 


2007-10-12 




















00053600031 


2007-10-24 




















00053600032 


2007-10-25 














SX3 






00053600033 


2009-10-16 


27 


1.0-8.0 


16392 


20 


5.0 


7.0 








00053600034 


2009-10-18 














XM1 


XMM 


MOS1 


0122700101 


2000-04-07 


29 


1.0-8.0 


77541 


200 


3.3 


5.0 


XM2 




MOS2 




2000-04-07 


29 


1.0-8.0 


76310 


200 


3.3 


5.0 


XPO 




pn 




2000-04-07 


24 


1.0-8.0 


172356 


200 


0.0 


2.75 



Notes. <n> Net exposure time after cleaning the events. (A) Energy range used in the spectral fitting. w Number of counts in the source extraction 
region in the fitting energy range. (rf) Number of counts per spectral bin for the fitting. {e) Inner and outer radii for the background extraction 
annulus. For ACIS-S3 data, we used a rectangular region (§ 14.2.11 . For EPIC-pn, we used blank sky data (§ 14. 2. 6b . 



a CdTe detector array of 128x128 pixels with a spectral reso- 
lution of ~9% at 100 keV. The IBIS has a coded mask, giving 
the instrument imaging capability of an unprecedented angular 
resolution of 12' in FWHM in this band. The full-coded field of 
view is ~9° square. 



dated systematic uncertainties, which could affect the analysis of 
bright sources. Due to the very long time span for the observa- 
tions of this source, we computed an averaged ARF for the total 
spectrum, which was weighted by the total exposure time in each 
calibration period. 



3.2.2. Data 

The IBIS data were obtained from a set of ~4000 IBIS Science 
Window data sets taken in a time period of 5 years from March 
2003 to March 2008. We checked the source light curve in the 
18-60 keV band and found that it is constant at a count rate of 
~0.9 s _1 . We thus stacked all the data sets. 

3.2.3. Processing 

The standard INTEGRAL Off-line Scientific Analysis (OSA) 
version 8.0 was used to obtain images and spectra. We used the 
standard IBIS response files released in the OSA 8.0 distribution. 
The IBIS ARFs are computed in each year to incorporate the up- 



3.3. RXTE PCA 

3.3.1. Instruments 

The PCA is a proportional counter array consisting of five pro- 
portional counter units (PCU1-5). The counters are filled with 
Xe gas and are sensitive at 2-60 keV. The instrument is unique 
in its high timing accuracy, a large effective area, and modest 
spectral resolution of ~1 keV at 6 keV. There are no focussing 
optics but the collimator limits the field of view to a ~1° square. 

3.3.2. Data 

We retrieved the archived data of G21.5-0.9 observed on 1997 
November 8 (table [TJ). We extracted the PCA spectrum from the 
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standard2 PCA data mode, which provides spectral information 
in 129 channels with a time resolution of 16 s. 

The electron flux stayed lower than 0.2 counts s _1 PCU 1 
and no PCU breakdown events occurred during the entire obser- 
vation. Data from all PCUs and all layers were combined into 
one spectrum. All five PCUs were operational during the en- 
tire observation, except for South Atlantic Anomaly (SAA) pas- 
sages. 

3.3.3. Processing 

We excluded events when the elevation angle from the Earth rim 
was smaller than 10° and during SAA passages and 30 minutes 
thereafter. As a result, we obtained a net exposure time of 21 ks. 

The background spectrum and the RMF were calculated us- 
ing the latest versions of pcabackest (version 3.8) and pcarmf 
(version 1 1 .7) tools respectively. The background spectrum was 
scaled down by 1.3% to match the Standard2 rate in channels 
110-128, corresponding to energies higher than 75 keV where 
no source counts are expected. 

3.4. Suzaku XIS and HXD-PIN 

3.4.1. Instruments 

The XIS is an X-ray imaging-spectrometer equipped with four 
X-ray CCDs sensitive in the 0.2-12 keV band. One is a back- 
side illuminated (XIS1) device and the others are front-side illu- 
minated (XIS0, 2, and 3) devices. The entire XIS2 and a part of 
the XIS0 are nonfunctional due to putative micro-meteorite hits, 
hence are not used in this paper. The four CCDs are located at 
the focal plane of four co-aligned X-ray telescopes with a HPD 
of ~2.'0. Each XIS sensor has 1024x1024 pixels and covers a 
17.'8xl7?8 view, which encompasses the entire G21.5-0.9 struc- 
ture and leaves ample room to construct background spectra. 

The PIN is a component of the HXD covering 10-70 keV. 
It is a non-imaging detector composed of 64 Si PIN diodes at 
the bottom of well-type collimators surrounded by GSO anti- 
coincidence scintillators. The effective area monotonically de- 
creases as the distance increases from the field center, with a full 
width at zero intensity (FWZI) view of ~70' square. 

3.4.2. Data 

We executed a Suzaku observation dedicated for the present 
study on 2009 October 10 for 40 ks using the XIS and the HXD. 
The target was placed at the XIS nominal position. The XIS was 
operated in the normal clocking mode with a frame time of 8 s. 

Within the FWZI square view of the PIN, two oth er sources 
were found in the INTEGRAL source catalogue ( Bir d et ail 
l20loh ; 4U 1835-11 and AX J1831.2-1008 respectively at 26'A 
and 42 .'7 from the field center. Because both of the unrelated 
sources have intensities less than 0.2% of G2 1.5-0. 9 in the 
ISGRI 20-60 keV band, we ignore the contribution of these 
sources in the PIN spectrum. 

3.4.3. Processing 

The data were processed through the pipeline processing version 
2.4 and were reduced using the HEASoft version 6.8-6.9. 

For the XIS, the RMFs and ARFs were generated using the 
xisrrafgen and xissimarfgen tools, respectively. The former 
takes into account the positional difference of the detector re- 
sponse. The latter is based on a ray-tracing simulation, which is 



designed to compensate for the lost photons outside of the pho- 
ton extraction region. The tool also takes into account the surface 
brightness distribution of incoming X-ray sources to calculate 
the effective area with the assumption that the spectrum is uni- 
form across the emitting region. We used the combined ACIS-S3 
image (figureQ} as the true spatial distribution. 

For the PIN, the background consists of astrophysical and in- 
strumental background. The former is dominated by the cosmic 
X-ray background (CXB). We used the instrumental background 
spectrum and the detector response files distributed by the instru- 
ment team. We simulated the CXB spectrum using the energy 
response and added it to the instrumental background to subtract 
the CXB contribution. G2 1.5-0. 9 is assumed to be point-like to 
generate the spatial response of the detector. 

3.5. Swift XRT 

3.5.1. Instruments 

The XRT comprises a Wolter-I telescope, which focuses X- 
rays onto an X-ray CCD device identical to t he ones flown 
on th e XMM-Newton EPIC-MOS instrument dBurrows et all 
l2004jh The CCD, which is responsive to 0.2-10 keV X-rays, has 
a dimension of 600x600 pixels, covering a 23'6x23.'6 field. The 
mirror has a HPD of ~18" and gives astrometric accuracies of a 
few arcseconds. 

As Swift's primary science goal is to rapidly respond to 
gamma-ray bursts (GRBs), the XRT was designed to operate 
autonomously, so that it could measure GRB light curves and 
spectra over seven orders of magnitude in flux. In order to mit- 
igate the effects of pile-up, the XRT automatically switches be- 
tween different CCD readout modes depending on the source 
brightness. Two frequently-used modes are: Windowed Timing 
(WT) mode, which provides ID spatial information in the cen- 
tral 7!8 of the CCD with a time resolution of 1.8 ms, and 
Photon Counting (PC) mode, which allows full 2D imaging- 
spectroscopy with a time resolution of 2.5 s. 

The XRT's thermoelectric cooler power supply system ap- 
parently failed shortly after launch, resulting in CCD operating 
temperatures ranging from -75 to -50 °C instead of the intended 
nominal -100 °C. In order to reduce the dark current and over- 
all noise level caused by the higher operating temperatures, the 
CCD substrate voltage (V ss ) was raised from to 6 V on 2007 
August 30. While this allowed the CCD to operate a few degrees 
warmer, it had the minor drawback of reducing the quantum ef- 
ficiency (QE) both at high energies (>5 keV) and just below the 
Si edge (1.5-1.84 keV). The V ss = 6 V calibration modifications 
required to correct the changes of QE have been completed for 
the WT mode but not for the PC mode, at the time of writing; 
the latter will be released in due course. 

3.5.2. Data 

G2 1.5-0. 9 is used as a routine calibration source for the XRT 
as its heavily absorbed spectrum is useful to verify the high en- 
ergy redistribution properties of the CCD. As the object is an 
extended source and it is difficult to subtract background ade- 
quately for the WT data, we present only results from PC mode 
observations. 

As Swift has a flexible observing schedule, often interrupted 
by GRBs, observations are divided into ObsIDs containing one 
or more snapshots, where each snapshot has a typical exposure 
of 1-2 ks. For a faint source like G2 1.5-0. 9, it is necessary to 
accumulate data from different ObsIDs for any given epoch. 
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We retrieved 18 observations from the UK Swift Science 
Data Centre. The observations were selected so that the source 
lay within 10' of the CCD bore-sight in order to avoid point- 
ings that have large off-axis angles to minimize the effect of vi- 
gnetting. The data were grouped into four epochs, split into two 
before (SXO and SX1) and after (SX2 and SX3) the V ss change 
(tableQ]). Typical exposures were 17-18 ks per epoch. 

3.5.3. Processing 

The data were processed with the Swift software version 3.5, 
which was made available with HEASoft version 6.8. 

For each epoch, a summed exposure map was computed, 
which accounts for the CCD bad columns and the telescope vi- 
gnetting. The exposure maps were then used to create extended 
source ARF files, by subdividing the source extraction region 
into boxes and calculating an average ARF over the regions, 
weighted by the relative exposure and counts in each box. For 
this analysis, the CALDB vOl 1 RMF and ARFs were used. 

3.6. XMM-Newton EPIC 
3.6.1. Instruments 

The EPIC is composed of two different types of X-ray CCD 
devices; two units of MOS-type (MOS1 and MOS2) covering 
0.15-12 keV and one unit of pn-type covering 0.15-15 keV. The 
three devices are located at the focal plane of three co-aligned 
telescopes with a HPD of ~15". 

The MOS 1 and MOS2 cameras consist of an array of seven 
front-illuminated CCD chips, each of which has 600x600 pixels 
covering a ~10f9xl0!9 square region. The pn-type has a 2x6 
back-illuminated CCD array, each of which has 200x64 pixels 
covering a ~13'6x4.'4 rectangular region. 
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Fig. 3. Layout of the field of view of the central chip of the two 
MOS arrays (solid square), the MOS and pn source extraction 
region (solid circle), and the MOS background extraction region 
(dashed annulus) overlaid on the pn image. The pn chip gaps and 
bad columns are apparent in the image. 



extraction region. However, we consider that this is negligible 
by taking a source extraction region large enough to encompass 
most photons to the level of 10~ 3 of the peak (figure[4]l. 



3.6.2. Data 

XMM-Newton observed G2 1.5-0. 9 several times within the 
field of view of the EPIC camera. We used the data with the tar- 
get placed at the center of the camera taken on 2000 April 7 (ta- 
ble[T}. The observation was performed with a medium-thickness 
filter and a full-frame mode for the CCD clocking. 

In the MOS arrays, the object is located at the center of the 
central CCD of both arrays. The entire spatial structure is con- 
tained in a chip. In the pn-array, G2 1.5-0. 9 is located at the bot- 
tom of CCD number 4 and the structure spreads across neigh- 
boring CCDs (figure©. 

3.6.3. Processing 

The data were reduced using th e XMM-Newton S cience 
Analysis System (SAS) version 10 dGabriel et al.ll2004l) using 
the most updated calibration available at the time of reducing 
the data (August 2010). Event lists were generated using the data 
reduction meta-tasks e[mp]proc with the default settings. The 
data were screened to remove intervals of high particle back- 
ground. 

Response files were generated using the rmf gen and arf gen 
tasks. We took into account the spatial distribution of the ex- 
tended X-ray emission by calculating a coarse detector map 
around the G2 1.5-0. 9 centroid, which was convolved with the 
spatially-dependent effective area to produce the appropriate 
transfer function for the extraction region. The SAS does not 
support the correction for the photons lost outside of the source 



4. ANALYSIS & RESULTS 

4.1. Overall procedure 

4.1.1. Source Extractions for Soft-band Instruments 

For the soft-band instruments, we extracted source events from a 
165" circle centered at the pulsar PSR J1833-1034 at (RA, Dec) 
= (18:33:33.57,-10:34:07.5) in the equinox J2000.0. The radius 
was chosen as a compromise (1) to be large enough to encom- 
pass a significant fraction of the entire structures of G2 1.5-0. 9 
as well as the negligible unrelated emission (figure|4]i and (2) to 
be small enough to be contained in the field of view of all the 
instruments. The Suzaku XIS with a larger telescope HPD than 
the others is the only exception; we chose a larger source ex- 
traction radius of 300", so that it contains a significant fraction 
of the source photons (figure 0). The background events were 
extracted from various regions optimized for each instrument, 
which is discussed individually in § 14.21 

4.1.2. Fitting Procedure 

We assembled the source and the background spectra, RMFs, 
and ARFs for all the instruments. Using the XSPEC package ver- 
sion 12.6, we fitted the background-subtracted spectra. We ap- 
plied the spectral model pegpwrlw atte nuated by the inters tellar 
photoelectric absorption model TBabs (Wil ms et alJ l2000). The 
former is the same with the conventional power-law model ex- 
cept that it uses the flux between the designated lower and upper 
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Fig. 4. Radial profile of background-unsubtracted 0.5-8 keV 
counts for the soft-band instruments. The counts are normalized 
to the value at the center. The radii of the source extraction aper- 
tures are shown with the dashed lines. 



energy as a free parameter for the intensity, not the normalized 
flux at 1 keV. This decouples the often entangled parameters of 
the index of power (F) and the normalization in power-law fit- 
ting. For the soft-band and hard-band instruments, the flux in 
the energy range of 2-8 k eV (Fx. so ft) and 15 -50 keV(Fx,hard) 
was evaluated, respectivelv. lVerner et al.l(ll996|) was used for the 
photoelectric absorption cross section, while lWilms et al.l (12000) 
was used for the metal abundance for the derivation of the H- 
equivalent extinction column density (Nh)- 

For the soft-band instruments, we used the common energy 
range of 1-8 keV for the fitting. All are Si-based CCDs, and thus 
have sensitivity in a similar energy range. For the hard-band in- 
struments, on the other hand, we used different ranges optimized 
for each instrument: 18-150 keV for INTEGRAL ISGRI, 5.0- 
30 keV for RXTE PC A, and 15-70 keV for Suzaku PIN. The 
inhomogeneity of the detecting technique did not allow us to use 
a common energy range. Because the power-law index of PWN 
spectra changes with the energy range, the choice of different fit- 
ting energy range might introduce an additional systematic un- 
certainty. 

Spectra were binned with the number of counts per bin spec- 
ified in tableQ] The model was fitted to minimize the^ 2 statistics 
and was considered statistically acceptable if the null hypothesis 
probability was larger than 5%. Statistical uncertainties for the 
best-fit values were derived as the 1 <x deviations. 



4.2. Individual Fitting and Systematic Uncertainties 

We now fit the spectra individually for each participating instru- 
ments following the common procedure set up in § 14.11 We jus- 
tify the choice of our spectral model by showing that all the fit- 
ting yields statistically acceptable results. We also discuss sys- 
tematic uncertainties arising from the common procedure. The 
spectra obtained from the same instrument and the same re- 
sponse (CSO-6, SX0-1, and SX2-3) are fitted jointly after con- 
firming that individual spectra give consistent results. 



4.2.1. Chandra ACIS 

The Chandra ACIS-S3 data (CS0-CS6) were taken at different 
epochs. G2 1 .5-0.9 was placed at different detector positions dis- 
placed 1 .'17—1 '32 from the chip center. The background spectra 
were accumulated from different regions (figure|2]l. Statistically- 
acceptable best-fit models were obtained, which are shown in 
figure [5] (the CX0 fitting as a representative) and table [2] 
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Fig. 5. ACIS-S3 (CS0) spectrum and the best-fit model. The up- 
per panel shows the data (crosses) and the best-fit model (solid 
line), while the lower panel shows the residuals to the fit. The 
range used for the fitting is shown with grey shading. 

In order to evaluate the systematic uncertainty caused by the 
inhomogeneities in the data sets, we compared the best-fit pa- 
rameters among the different observations (figure [6]). For each 
parameter of A^h, T, and Fx.soft, the seven values were tested 
against the null hypothesis that they deviate from a constant 
value. For all parameters, the hypothesis was rejected, imply- 
ing that there is no statistical reason to claim for any difference 
among these data sets. We thus combined all the spectra and con- 
ducted a joint fitting to obtain a statistically-acceptable best-fit 
model. The result is appended in table|2]as CSO-6. 

Due to the superb imaging capability of ACIS, the full-frame 
data suffer slight pile-up at the peak of the surface brightness dis- 
tribution of the object. We consider that this effect is negligible 
for our cross-calibration comparison based on the following as- 
sessment. We used the sub-array data and compared the results 
of the full-frame and sub-array fitting at different outer source 
extraction radii (r out ) within the field of view of the sub-array 
data (figure |7]i. A shorter frame time for the sub-array data re- 
duces the frequency of pile-up. In figure [7] we see two trends 
as r out increases. The first is the monotonic changes of the spec- 
tral shape (A^h and F) both for the full-frame and sub-array data, 
which stems from the spatial changes of the spectrum intrinsic 
to G2 1.5-0. 9. The second is the decreasing differences between 
the full-frame and sub-array results, which are due to the alle- 
viation of the pile-up effect. Beyond 20", the results of the two 
data sets are stable with a full-frame to sub-array ratio of 1 ex- 
cept for a 4% difference in Fx, so ft- We consider that the Fx, S oft 
difference is a separate issue with an unidentified origin. In fact, 
we found a similar discrepancy in a comparison of an annular 
region, in which the inner 30" circle including the bright core 
was removed to avoid pile-up. The Fx, S oft difference between the 
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Table 2. Best-fit parameters and uncertainties in the spectral fitting. 



Label 




r 


^X.soft 


-Fx.hard'' 


Red-^ 2/ 




(10 22 cm- 2 ) 




(1Q- 11 ergs 


-' cm- 2 ) 


/d.o.f. 


Chandra ACIS-S3 


CSO 


2.99 (2.93-3.04) 


1.83 (1.80-1.86) 


6.10(6.05-6.16) 




0.93/ 302 


CS1 


3.07 (3.01-3.12) 


1.85 (1.83-1.88) 


6.09 (6.04-6.13) 




0.90/ 326 


CS2 


3.04 (2.98-3.09) 


1.82 (1.79-1.84) 


6.06 (6.01-6.11) 




1.04/325 


CS3 


3.11 (3.05-3.16) 


1.84(1.81-1.87) 


6.04 (5.99-6.09) 




0.89/ 327 


CS4 


3.16(3.11-3.22) 


1.88 (1.85-1.91) 


6.10(6.05-6.15) 




1.03/330 


CS5 


3.00 (2.95-3.06) 


1.81 (1.78-1.84) 


6.01 (5.97-6.06) 




1.06/327 


CS6 


3.14 (3.08-3.20) 


1.88 (1.85-1.91) 


6.03 (5.98-6.08) 




1.07/326 


CSO-6 


3.07 (3.05-3.09) 


1.84(1.83-1.85) 


6.06 (6.04-6.08) 




0.99/2281 


INTEGRAL IBIS-ISGRI 


ISO 


2.99 


2.09 (2.02-2.16) 




4.17 (4.03-4.31) 


1.43/9 


RXTE PCA 


RPO 


2.99 


2.05 (2.04-2.07) 




5.54 (5.45-5.64) 


0.61/51 


Suzaku XIS 


SIO 


3.07 (3.03-3.11) 


1.91 (1.89-1.92) 


5.62 (5.59-5.65) 




1.04/616 


SI1 


3.18(3.14-3.22) 


1.92 (1.90-1.93) 


5.73 (5.70-5.76) 




1.06/632 


SI3 


3.07 (3.03-3.11) 


1.90(1.88-1.92) 


5.66 (5.63-5.69) 




0.98/ 624 


Suzaku HXD-PIN 


SPO 


2.99 


2.28 (2.14-2.42) 




6.10 (5.79-6.42) 


1.40/ 12 


Swift XRT 


SXO 


2.97 (2.88-3.07) 


1.77 (1.73-1.81) 


5.79 (5.72-5.87) 




0.99/421 


SX1 


2.90 (2.83-2.98) 


1.77 (1.74-1.81) 


5.48 (5.42-5.54) 




1.03/479 


SX2 


3.05 (2.98-3.13) 


1.90(1.87-1.94) 


5.46 (5.40-5.51) 




1.07/488 


SX3 


3.16(3.08-3.25) 


1.93 (1.89-1.96) 


5.46 (5.40-5.52) 




1.14/478 


SX0+1 


2.93 (2.87-2.99) 


1.77 (1.75-1.80) 


5.61 (5.56-5.65) 




1.02/903 


SX2+3 


3.10(3.05-3.16) 


1.91 (1.89-1.94) 


5.46 (5.41-5.50) 




1.11/ 969 


XMM-Newton EPIC 


EMI 


2.90 (2.87-2.94) 


1.78 (1.77-1.80) 


5.47 (5.44-5.50) 




1.05/264 


EM2 


2.94 (2.91-2.98) 


1.84(1.83-1.86) 


5.37 (5.34-5.40) 




1.11/ 260 


EPO 


2.74 (2.72-2.77) 


1.76(1.75-1.77) 


5.07 (5.06-5.09) 




1.13/623 


Joint fitting 


Soft 


2.99 (2.98-3.00) 


1.84(1.84-1.85) 


5.69 




1.07/6217 


Hard 


2.99 


2.05 (2.01-2.09) 




4.87 


0.96/ 66 



Notes. <D) The parentheses indicate the 1 tr statistical uncertainty of the best-fit value. The values without the uncertainty indicate that the values 
were fixed in the fitting. <6) H-equivalent column density of the interstellar extinction. w The power-law index. (rf) Flux in the 2.0-8.0 keV range 
for the soft-band instruments. <c) Flux in the 15-50 keV range for the hard-band instruments. The goodness of fit with the reduced^ 2 (Red-^ 2 ) 
and the degree of freedom (d.o.f.). 
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Fig. 6. Comparison of best-fit values and their 1 <x statistical un- 
certainty for the spectral parameters by the full-frame data of 
ACIS-S3. The solid line and the grey shaded regions respectively 
indicate the best-fit and 1 cr statistical uncertainty range of joint 
fitting. 



full-frame and sub-array data was again found at a ~3.5% level, 
which exceeds the statistical uncertainty of ~2.0%. We thus con- 
clude that the major cause of the Fx,soft discrepancy is not pile- 
up. 



4.2.2. INTEGRAL IBIS-ISGRI 

The IBIS-ISGRI detected emission from the object at a signif- 
icance level of ~35cr. The spectrum consists of 12 channels 
spaced almost equally in a logarithmic scale (figure|S). The inter- 
stellar extinction cannot be constrained in the IBIS-ISGRI data 
alone, so we fixed the value to 2.99xl0 22 cirr 2 , which was ob- 
tained in the joint fitting of the soft-band instruments (§ 14. 3. U . 
A statistically-acceptable best-fit model was obtained, which is 
shown in figure [8] and table [2] 

In order to compare with the previous work bv lde Rosa et al.l 
(2009) for the same object, we fitted the spectrum with the same 
model (a single power-law) used in their work. The derived index 
and flux in the 20-100 keV range is 2.0 (1 .9-2.1) and 5.3 (5.1- 
5.5) xl0~ n erg s _1 cirT 2 , whereas those bv lde Rosa et alJ (T2009) 
are 2.2 + 0.1 and 5.2 xlO -11 erg s -1 cirr 2 . They are consistent 
with each other within the statistical uncertainty. 
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2.99xl0 22 cirT 2 similarly to the IBIS-ISGRI fitting. All param- 
eters in the model for the unrelated emission were common be- 
tween the two spectra, except for the flux rescaling factor for the 
G2 1.5-0. 9 spectrum with respect to the control field spectrum. 

As a result, statistically-acceptable best-fit models were ob- 
tained (figure [9]). The best-fit parameters for the G2 1.5-0. 9 
model are summarized in table [2] Those for the unrelated emis- 
sion model are; T =2.39 (2.25-2.51), F XMrd =1.50 (1.24-1.83) 



x 10 11 erg s 1 cm 2 , E gml 
(0.34-0.45) keV, and JV gau 



=6.56 (6.53-6.59) keV, <x g 
=3.85 (3.57-4.18) xlO 4 s _1 cm 



0.40 

2 



where E gau , o" gau , and N g . du are the center energy, width, and flux 
of the Gaussian model. In the 15-50 keV range, it is estimated 
that the contribution of the unrelated emission accounts for 21% 
of the total emission in the G2 1.5-0. 9 spectrum. Note that the un- 
certainties in the control-field flux was not propagated to derive 
the numbers in table [2] 



Fig. 7. Comparison of best-fit values and their 1 <x statistical 
uncertainty for the spectral parameters by joint fitting of the 
full-frame (black) and the sub-array (red) data at different outer 
source extraction radii. The inner radius was kept 0. The ratio of 
the two (green) are also shown for the right horizontal axis. 
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Fig. 9. PCA spectrum and the best-fit model. In the upper panel, 
the best-fit model (thick brown solid line) is decomposed into the 
G21.5-0.9 (brown dotted line) and the unrelated (orange dotted 
200 line) emission components. The other symbols follow figure[5] 



Fig. 8. ISGRI spectrum and the best-fit model. The symbols fol- 
low figure [5] 

4.2.3. RXTE PCA 

The non-imaging nature and the large field of view of the PCA 
lead to a considerable contamination by unrelated emission ubiq- 
uitous in the direction of G2 1.5-0. 9. Two major sources of con- 
tamination are: (1) the Galactic ridge X-ray emission, which is 
evident from an excess emission at 6-7 keV (figure [9} and (2) 
the CXB emission. 

In order to subtract the contribution of the unrelated emis- 
sion, we used the data (ObsID = 20266-01-35-06S) obtained in 
a control field devoid of bright sources at the Galactic coordinate 
of(/, b) = (21.999, -0.005), which is 1°.01 away from G2 1.5-0.9 
at (21.501, -0.885). We fitted the two spectra at the control field 
and at G2 1.5-0. 9 jointly. For the former, a model was applied to 
represent the unrelated emission, which is a power-law contin- 
uum plus a Gaussian line for the blended Fe lines at 6-7 keV at- 
tenuated by a fixed interstellar extinction of Nu =3x10 22 cm" 2 . 
For the latter, the model representing the G2 1.5-0. 9 spectrum 
was added. The interstellar extinction to G2 1.5-0. 9 was fixed to 



4.2.4. Suzaku XIS and HXD-PIN 

The background events were extracted from an annulus of inner 
and outer radii of 5'.0 and 7'0, which is concentric to the source 
extraction region. The calibration of the XIS is uncertain around 
the Si K edge. We excluded the events in the 1.8-2.0 keV range 
in the fitting. We performed independent fitting for each of the 
three XIS and a PIN spectra. A statistically-acceptable best-fit 
models were obtained, which is shown in figure[T0]and table|2] 

Unlike the other soft-band instruments, the ARF generator of 
the XIS compensates for lost photos outside of the source extrac- 
tion region. This is a source of systematic uncertainty. In order 
to evaluate this, we fitted the spectra constructed from different 
outer source extraction radii (r out ). The inner radius was fixed to 
0'. The ARFs were generated by assuming a point-like source 
or a surface brightness distribution of figure Q] Figure QT| sum- 
marizes the result, in which we see the best-fit spectral param- 
eters converge as the increasing r out for both assumptions. For 
'"out =5', which we use to compare with the other instruments, 
we can safely assume that the uncertainty due to the compensa- 
tion of lost photons is negligible. 
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Fig. 10. XIS and PIN spectra and the best-fit model. The symbols 
follow figure [5] 
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Fig. 11. Best-fit values and their 1 cr statistical uncertainty for the 
spectral parameters with different outer source extraction radii. 
Black symbols indicate that the ARF was calculated assuming 
the source brightness distribution of the ACIS-S3 image. Red 
symbols indicate that a point-like source was assumed. Here, we 
use the result of joint fitting of the three XIS devices upon con- 
firmation that the result of individual fitting gives the same trend. 



4.2.5. Swift XRT 

Figure[T2] shows the summed XRT image from which the radial 
profile (figure [4]i was calculated. Overplotted on the image are 
the source (solid) and background (dashed) extraction regions. 
The background was chosen from an annuls and some additional 
circles to increase the area from the region with a maximum and 
uniform exposure. 

For each of the four epochs of data (table [TJ, we performed 
individual fitting. We further grouped all the data sets before and 
after the V ss change for joint fitting. Statistically-acceptable best- 
fit models were obtained, which are shown in figure [5] (the SX0 
fitting as a representative) and table [2] 

There is a clear difference in the spectral parameters obtained 
from data taken before and after the V ss change (figure[T4"l). with 
the latter showing a slightly steeper but more absorbed spectrum 
than the former. Within each group, Nh and T are consistent with 



28:00 



40:00 i 




18:34:00.0 50.0 40.0 30.0 20.0 10.0 

RA (J2000.0) 

Fig. 12. XRT image of G2 1.5-0. 9. The source (solid) and back- 
ground (dashed) regions used for the spectral analysis are shown. 




1.0 _ 2.0 

Energy (keV) 



Fig. 13. XRT (SX0) spectrum and the best-fit model. The sym- 
bols follow figure[5] 



each other. The discrepancy in the parameters is caused by the 
slight QE change incurred after the V ss was altered. The cali- 
bration of this effect is ongoing. We hereafter use the result ob- 
tained for the data before the change (SX0 and SX1). We see a 
notable difference in fx, S oft between SX0 and SX1 . We could not 
identify the cause for the difference and ignore it as a statistical 
scatter. 



4.2.6. XMM-Newton EPIC 

In order to ensure that the fluxes measured by the three EPIC 
cameras correspond to the same region in the sky, as well as to 
provide a correction factor for the lost sky area due to bad pix- 
els/columns and chip gaps, we applied a common detector mask 
during the accumulation of the spectral products ("masked spec- 
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£ SX1 SX2 SX3 SX4 0.5 1.0 2.0 5.0 

Data group Energy (keV) 



Fig. 14. Comparison of best-fit values and their 1 cr statistical 
uncertainty for the spectral parameters by XRT. The solid line 
and the grey shaded regions respectively indicate the best-fit and 
1 cr statistical uncertainty range of joint fitting. 



tra" hereafter). This mask was generated by multiplying the de- 
tector masks produced by the SAS task emask on the individual 
EPIC camerafl For the correction factor to compensate for the 
masked area, we extracted a MOS1 spectrum from the full 165" 
extraction circle ("unmasked spectrum"). We chose the MOS1 
detector here, because it was only marginally affected by CCD 
artefacts in our data. We calculated the flux ratio between the un- 
masked and the masked spectra, which was 1.091 in the 2-8 keV 
energy range. Finally, we multiplied the ratio to the flux obtained 
by the "masked" spectra to derive the flux in the full 165" circle. 

The background spectra for the two MOS cameras were ex- 
tracted from an annulus of the inner and outer radii of 33 and 
5'0, which is concentric to the source extraction region. The 
background region is also contained in the central CCD chip 
with the source region (figure [3]). 

Unlike the MOS cameras, the background extraction for the 
pn camera is not straight-forward. The source spreads across 
multiple CCDs (figure [3} and the background level is known to 
have a measurable position dependence, in particular, along the 
readout di rection. We therefor e extracted a spectrum from blank 
sky data (ICarter & R ead 2007) using the same extraction region 
with the source and used it as the background spectrum. Out-of- 
time events were removed in the pn data. 

We fitted individual spectra of the three EPIC devices. A 
statistically-acceptable best-fit model was obtained, which is 
shown in figure Q3] and table [2] 

In the pn result, a deviation is seen below 1.5 keV, which is 
not observed in any other instruments at a comparable level. The 
origin of this feature is currently unknown. Possible explanations 
are residual uncertainties in the background subtraction (the pn 
background has a complex spatial structure whose modelling for 
extended sources close to the bore-sight position is non-trivial) 
and residual inaccuracies in the calibration of the redistribution. 



2 SAS users should be aware that this procedure does not work with 
SAS versions 10 and earlier. Due to a formal error in the header of 
the mask file, the mask description is not correctly propagated in the 
spectral file header, and arfgen cannot properly apply the bad area 
correction on "masked" spectra. We modified the header of the mask file 
manually in order to cope with this problem. A solution to this problem 
is expected to be available in SAS versions later then 10.0. 



Fig. 15. EPIC spectra and the best-fit model. The symbols follow 
figure [5] 



Notwithstanding the origin of this feature, we assessed the effect 
of the deviation to the best-fit parameter values by altering the 
lower energy (.Emm) in the fitting (figure [T6l>. The upper energy 
(£max) was fixed to 8.0 keV. Toward lower energies in E m m, the 
best-fit A^h and F values change monotonically. The global trend, 
however, is common among MOS1, MOS2, and pn, implying 
that the <1.5 keV deviation only seen in the pn does not affect 
the pn fitting results. This is conceivable given the fact that the pn 
spectrum has dominant counts in 1 .5—E max keV band. We thus 
ignore the effect of the <1.5 keV deviation hereafter. 



W E 4.5 
3 4.0 

CM 

w o 3.5 

~3.0 

2 2.5 
2.1 

2.0 

U 1.9 

_ 1-8 

°E 1.7 

5.8 

I 5-6 

- 5.4 

° 5.2 

= 5.0 



pn 



m m w 



8 8 8 



nut 

-I 1 1 1 H 



„ CD * <t <j> * 



CD CD CD CD 



* ♦ + M M M 



_, I ._ 



0.5 



1.0 



1.5 



2.0 



2.5 



: lower energy band in fitting (keV) 



Fig. 16. Best-fit values and their 1 cr statistical uncertainty for 
the spectral parameters with different lower energy band used in 
the fitting. 



4.3. Joint Fitting of All Data 

We assembled all the spectra and responses and attempted a joint 
fitting separately for the soft-band and hard-band instruments. 
The procedures follow the individual fitting in § 14.21 The seven 
ACIS-S3 spectra were treated as one. The two XRT data (SX2 
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and SX3) were discarded and the remaining two (SXO and SX1) 
were used as one. 



4.3.1. Soft-band Instruments 

In the soft-band joint fitting, we first fitted the spectra with all 
the spectral parameters (A^h, T, Fx.sott) tied among all the instru- 
ments. The best-fit model was statistically rejected due mainly 
to the inconsistent normalizations among the instruments. The 
statistical uncertainties were not derived as the goodness of fit 
was very low. 

We thus introduced an additional normalization rescaling pa- 
rameter for all the instruments with respect to a fixed Fx.soft 
value of 5.69X1CT 11 erg s _1 crrT 2 , which is the best-fit value in 
the joint fitting above. Note that this value is not the most likely 
value for the absolute flux. The derived value can be easily bi- 
ased by the instruments yielding a larger number of counts than 
the others regardless of whether they are well-calibrated or not. 
Nevertheless, it is a value that best represents our datasets. We 
obtained a result, which was still statistically unacceptable but 
was much improved in the goodness of fit (table The com- 
parison of the F x.soft is shown in figure [171 
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Fig. 18. Scatter plot of the best-fit F versus A^h values and their 
1 cr statistical uncertainty of all soft-band instruments. 
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Fig. 17. Comparison of the best-fit Fx, so ft or Fx.hard values and 
their 1 cr uncertainty in the joint fitting with the other parameters 
(T and Nh) tied among all the soft-band or hard instruments. 
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4.3.2. Hard-band Instruments 

We took the same approach in the joint fitting of the hard- 
band instruments except that the Nh value was fixed to 
2.99xl0 22 crrT 2 , which is the result of the soft-band joint fit- 
ting. First, we fitted all the spectral parameters (T, Fx.hard) tied 
among all the instruments, which was in vain for the same rea- 
son with the soft-band joint fitting. We thus introduced the nor- 
malization rescaling parameter for each spectrum and obtained a 
statistically-acceptable best-fit model (table[2]i. The comparison 
of the Fx.hard is shown in figure [T71 



Fig. 19. Scatter plot of the best-fit A% versus Fx, S oft values and 
their 1 cr statistical uncertainty of all soft-band instruments. 

5. DISCUSSION 

5. 1 . Comparison of Fitting Results 

We now compare all the results. The values in tableware plotted 
in four scatter plots (figure [18] [19] [20] and [21]) for every combi- 
nation of two spectral parameters. 
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Fig. 20. Scatter plot of the best-fit F versus Fx.soft values and 
their 1 cr statistical uncertainty of all soft-band instruments. 
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Fig. 21. Scatter plot of the best-fit T versus Fx.hard values and 
their 1 <x statistical uncertainty of all hard-band instruments. 

5.1.1. Identifying Systematic Differences 

For identifying significant systematic uncertainties, we assume 
that the instrument A and B respectively have the best-fit value 
and 1 cr statistical uncertainty for a certain spectral parameter 



x as x (A) + Ajc (A) and 



,(B) 



Ajc (b \ The statistical uncertainty is 



nearly symmetric in the upper and lower directions in all param- 
eters (table 0, so we took the mean of the two for Ax (A) and 



Ax (B) here. If the modulus of the logarithm of the ratio exceeds 
three times the convolved statistical uncertainty as 



In 



> 3 



//Ax( A >\ 2 /Ax< B >\ 2 



(1) 



we consider that the two instruments A and B have a significant 
systematic difference in x. 

In the soft band, we have eight instruments (ACIS-S3, XIS0, 
XIS1, XIS3, XRT, MOS1, MOS2, and pn) and three parameters 
(A^h, T, and Fx.soft)- Figure[T8]shows that A^h and F are coupled, 
so we only consider T. For each parameter, we have 28 combi- 
nations for a pair of instruments. In the hard band, we have three 
instruments (PIN, PCA, and ISGRI) and two parameters (r and 
Fx.hard)- For each parameter, we have 3 combinations. A total of 
62 combinations were examined. The identified systematic dif- 
ferences are shown in bold faces in tables [3j SJ E] and[6] 

5.1.2. Soft-band Comparison 

In the flux comparison, 20 out of 28 combinations show a sig- 
nificant difference (table [3j. The ACIS-S3 flux is larger than all 
the others. The pn flux is smaller than all the others. The three 
XIS instruments (XIS0, XIS1, and XIS3) are consistent with 
each other, while the two MOS cameras are consistent with each 
other. XRT is consistent with the three XIS devices and MOS 1 . 
The largest difference, between ACIS-S3 and pn, is 20%. 

In the index comparison, 16 out of 28 combinations show a 
significant difference (table |4|. The three XIS instruments and 
the two MOS cameras are consistent with each other also in T. 
The largest difference, between XIS1 and pn, is 9%. 



5.1.3. Hard-band Comparison 

In the flux comparison, 2 out of 3 combinations in fx.hard are 
significantly different (table [5j. ISGRI is smaller than the other 
two. The largest difference between ISGRI and PIN is as large 
as 46%. In the index comparison, all three instruments are con- 
sistent with each other (table[6j. 



5.2. Comparison with Other Studies 

5.2.1. Internal Cross-calibration Studies 

We compared our results with the cross-calibration results for 
the instruments onboard the same satellite (i.e., XIS0, XIS1, and 
XIS3 for Suzaku and EPIC-MOS1, EPIC-MOS2, and EPIC-pn 
for the XMM-Newton), referring to the technical notes issued 
by their calibration teams. Both observatories use these instru- 
ments simultaneously, so we can mitigate the effects caused by 
the changes of the instrumental performance in time in these 

comparisons. 

In the Suzaku XIS study dlshida et al]l2007l) . the Crab neb- 
ula was used to compare the power-law index and the flux in 
1-10 keV at the XIS nominal position. Their results are consis- 
tent with the present study. XIS0, 1, 3 are consistent with each 
other to better than 3% both in the flux and power-law index. The 
back-side illuminated device (XIS1) is slightly larger in flux and 
steeper in the power-law index than the front-side illuminated 

devices (XIS0 and XIS3). 

In the XMM-Newton EPIC study dGabriel et alJ 120081) . 
G21.5-0.9 was also used to compare the power-law index and 
the flux in 2-10 keV band using the same data set. Their results 
shows that the flux increases in the sequence of MOS2, MOS1, 
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Table 3. Relative flux among soft-band instruments." 





A 




ACIS-S3 


XISO 


XIS1 


XIS3 


XRT 


MOS1 


MOS2 


pn 




ACIS-S3 


0.00 


7.54 ± 0.63 


5.60 ± 0.62 


6.83 ± 0.62 


7.72 ± 0.87 


10.24 ± 0.64 


12.09 ± 0.65 


17.84 ± 0.44 




XISO 


-7.54 ± 0.63 


0.00 


-1.94 ± 0.75 


-0.71 +0.75 


0.18 ±0.96 


2.71 ± 0.77 


4.55 ± 0.77 


10.30 ± 0.61 




XIS1 


-5.60 ± 0.62 


1.94 ± 0.75 


0.00 


1.23 ± 0.74 


2.12 ±0.96 


4.64 ± 0.76 


6.49 ± 0.77 


12.24 ± 0.60 




XIS3 


-6.83 ± 0.62 


0.71 ± 0.75 


-1.23 ± 0.74 


0.00 


0.89 ± 0.96 


3.41 ± 0.76 


5.26 ± 0.77 


11.01 ± 0.61 


B 


XRT 


-7.72 ± 0.87 


-0.18 ± 0.96 


-2.12 ±0.96 


-0.89 ± 0.96 


0.00 


2.53 ± 0.97 


4.37 ± 0.98 


10.12 ± 0.85 




MOS1 


-10.24 ± 0.64 


-2.71 ± 0.77 


-4.64 ± 0.76 


-3.41 ± 0.76 


-2.53 ± 0.97 


0.00 


1.85 ± 0.78 


7.59 ± 0.62 




MOS2 


-12.09 ± 0.65 


-4.55 ± 0.77 


-6.49 ± 0.77 


-5.26 ± 0.77 


-4.37 ± 0.98 


-1.85 ± 0.78 


0.00 


5.75 ± 0.63 




pn 


-17.84 ± 0.44 


-10.30 ± 0.61 


-12.24 ± 0.60 


-11.01 ± 0.61 


-10.12 ± 0.85 


-7.59 ± 0.62 


-5.75 ± 0.63 


0.00 



Notes. <n) The logarithmic flux ratio 100 x ln^F^V^jf*) between the soft-band instruments A and B. The flux is measured in the 2.0-8.0 keV 

band. The range indicates the convolved statistical uncertainty derived as lOOx aN "psrj + y^W J ■ Both the ratio and its deviation are multiplied 

by 100 to save space. The uncertainties AF^' and AF^ 5 ' are the mean of the 1 a statistical uncertainties in the upper and lower bound directions 
(table |2} . The bold face indicates that the difference is larger than 3 times the convolved statistical uncertainty. 



Table 4. Relative index of power among soft-band instruments. 13 





A 


ACIS-S3 


XISO 


XIS1 


XIS3 


XRT 


MOS1 


MOS2 


pn 




ACIS-S3 


0.00 


-3.73 ± 0.96 


-4.26 ± 0.95 


-3.21 ± 1.18 


3.88 ± 1.51 


3.32 ± 1.00 


0.00 ± 0.98 


4.45 ± 0.79 




XISO 


3.73 ± 0.96 


0.00 


-0.52 ±1.11 


0.52 ± 1.31 


7.61 ± 1.62 


7.05 ± 1.15 


3.73 ± 1.13 


8.18 ± 0.97 




XIS1 


4.26 ± 0.95 


0.52 ± 1.11 


0.00 


1.05 ± 1.31 


8.13 ± 1.61 


7.57 ± 1.15 


4.26 ± 1.13 


8.70 ± 0.97 




XIS3 


3.21 ± 1.18 


-0.52 ± 1.31 


-1.05 ± 1.31 


0.00 


7.09 ± 1.76 


6.52 ± 1.35 


3.21 ± 1.33 


7.65 ± 1.20 


B 


XRT 


-3.88 ± 1.51 


-7.61 ± 1.62 


-8.13 ± 1.61 


-7.09 ± 1.76 


0.00 


-0.56 ± 1.64 


-3.88 ± 1.63 


0.57 ± 1.52 




MOS1 


-3.32 ± 1.00 


-7.05 ± 1.15 


-7.57 ± 1.15 


-6.52 ± 1.35 


0.56 ±1.64 


0.00 


-3.32 ± 1.17 


1.13 ± 1.02 




MOS2 


0.00 ± 0.98 


-3.73 ± 1.13 


-4.26 ± 1.13 


-3.21 ± 1.33 


3.88 ± 1.63 


3.32 ± 1.17 


0.00 


4.45 ± 0.99 




pn 


-4.45 ± 0.79 


-8.18 ± 0.97 


-8.70 ± 0.97 


-7.65 ± 1.20 


-0.57 ± 1.52 


-1.13 ± 1.02 


-4.45 ± 0.99 


0.00 



Notes. (a) The logarithmic index of power ratio 100 x ln(r <B) /T (A) ) between the soft-band instruments A and B. The notation follows tablefj] 



Table 5. Relative flux among hard-band instruments." 





A 




ISGRI 


PCA 


PIN 




ISGRI 


0.00 


-28.41 ± 3.77 


-38.04 ± 6.16 


B 


PCA 


28.41 ± 3.77 


0.00 


-9.63 ± 5.44 




PIN 


38.04 ± 6.16 


9.63 ± 5.44 


0.00 



Notes. <n) The logarithmic flux ratio 100 x lnfF^ /F x J between the hard-band instruments A and B. The flux is measured in the 15-50.0 keV 
band. The notation follows tablef3] 



and pn, which is inconsistent with ours. This is due to the caveat 
of the bad pixel/column treatment in the older SAS versions de- 
scribed in § 14.2.61 In fact, w ithout considering t his caveat, we 
obtained the same result with lGabriel et alJ (120081) . 

In another XMM-Newton EPIC study dStuhlinger et all 
2008), several extra-galactic point-like sources are used to com- 
pare the flux at a half dozen energy bands. Point sources are less 
affected by the SAS caveat than extended sources like G21.5- 
0.9. In the two energy bands closest to ours (1.5^.0 and 4.0- 
10.0 keV), the flux is larger in MOS 1 and M02 than pn by -5%, 
which is the same trend with our result. 



5.2.2. External Cross-calibration Studies 

We also compared our results with other cross-calibration stud- 
ies across multiple missions. The number of such studies is ac- 
tually small; in fact, we could only find several studies compar- 
ing Chandra and XMM-Newton often with outdated software 



versio ns. We refer to another IACHEC study iNevalainen et alJ 
(1201 Oh . in which 1 1 relaxed clusters of galaxies are used. 

The 2-7 keV band spectra were fitted with an optically-thin 
thermal plasma model to derive the plasma temperature and the 
flux. The ACIS-I, ACIS-S, and EPIC instruments were used. 
ACIS-I and ACIS-S were confirmed to be consistent with each 
other and were treated as one. 

In the flux comparison, ACIS was significantly larger than 
the three EPIC cameras by 5-10%. Among the three EPIC in- 
struments, pn is lower than the other two by ~5%. Both discrep- 
ancies are consistent with our result. 

In the temperature comparison, ACIS and pn were consis- 
tent with each other. Among the EPIC instruments, the pn and 
the MOS (MOS 1 and MOS2 were combined) were also consis- 
tent with each other. In order to compare to the cluster temper- 
ature result, we fitted the G2 1.5-0. 9 spectra in the same energy 
band with a power-law model attenuated by a fixed A^h value of 
2.99xl0 22 ctrr 2 and substituted the softness of the index with 
the temperature. We found that the consistency in the spectral 
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Table 6. Relative index of power among hard-band instruments." 





A 


ISGRI 


PCA 


PIN 




ISGRI 


0.00 


1.93 ± 3.43 


-8.70 ± 6.99 


B 


PCA 


-1.93 ±3.43 


0.00 


-10.63 ±6.18 




PIN 


8.70 ± 6.99 


10.63 ± 6.18 


0.00 



Notes. (fl) The logarithmic index of power ratio 100 x ln(r <B) /T (A) j between the hard-band instruments A and B. The notation follows table[3] 



softness among the ACIS-S3 and the three EPIC cameras is also 
found in G21. 5-0.9. 

6. CONCLUSION 

We used the pulsar wind nebula G2 1.5-0. 9 to conduct a cross- 
calibration study of instruments onboard the currently working 
X-ray astronomy missions. The archived as well as the origi- 
nal data were accumulated for Chandra ACIS-S3, INTEGRAL 
IBIS-ISGRI, RXTE PCA, Suzaku XIS and HXD-PIN, Swift 
XRT, and XMM-Newton EPIC-MOS and EPIC-pn. 

We conducted a coherent spectral fitting for all the data with 
the methods as common as possible. We tabulated and plotted 
the results to be a useful reference. We compared the results in all 
combinations of two different instruments in the spectral param- 
eters and identified many systematic differences unattributable 
to the statistical uncertainties alone. Most differences are consis- 
tent with those reported in previous studies of a similar purpose. 

Many factors can cause the systematic differences among 
different instruments, including (1) inaccurate calibration of ei- 
ther or both of the two compared instruments, (2) the limita- 
tions of G2 1.5-0. 9 as a calibration source, and (3) inappropri- 
ate choice of the source extraction aperture, the spectral models, 
background, energy ranges, etc. These are entangled with each 
other, and our aim is to discriminate the causes belonging to (1). 
We addressed some of the concerns belonging to (3) in § !4.2l and 
showed that their effects are minor. The limitations of G2 1.5-0. 9 
as a calibration source (§[2]) are unavoidable. After all, there is no 
perfect celestial calibration source. Therefore, the only practical 
approach is to conduct similar comprehensive comparison stud- 
ies using a variety of celestial sources and to identify systematic 
differences commonly seen in various studies. More IACHEC 
papers will follow, which will serve for this purpose. 

The systematic differences should not be left as they are. We 
feel that the current level of differences among various instru- 
ments (or even within the same instrument) can be narrowed to 
a greater extent. For example, we described in this study that (1) 
the Chandra ACIS-S3 full-frame and the sub-array data appear 
to have inconsistent normalization, (2) the Swift XRT results re- 
main inconsistent before and after the substrate voltage changes 
or even at two epochs before the change for flux, (3) Suzaku 
XIS is not well calibrated at the energies around the Si I K edge, 
(4) the EPIC-pn shows a deviation in the low energy tail of the 
spectrum, and (5) the three EPIC instruments are inconsistent 
with each other. 

Scientific measurements can be no better than the calibration 
accuracies. Continuous calibration efforts are requisites to keep 
the instruments reliable and to maximize the scientific output. 
We hope that the present study will help such efforts, which the 
entire community relies upon. 

Acknowledgements. We thank the IACHEC members for sharing their results 
prior to publication. We acknowledge Dai Takei and Kei Saitou for their help in 
Suzaku data reduction. APB and AMR acknowledges support from the STFC. 



References 

Asaoka, I. & Koyama, K. 1990, PASJ, 42, 625 

Bandiera, R. & Bocchino, F. 2004, Advances in Space Research, 33, 398 
Becker, R. H. & Szymkowiak, A. E. 1981, ApJ, 248, L23 
Beuermann, K. 2006, A&A, 460, 783 
Bietenholz, M. F. & Bartel, N. 2008, MNRAS, 386, 1411 
Bird, A. J., Bazzano, A., Bassani, L., et al. 2010, ApJS, 186, 1 
Bocchino, F., van der Swaluw, E., Chevalier, R., & Bandiera, R. 2005, A&A, 
442, 539 

Bradt, H. V, Rothschild, R. E., & Swank, J. H. 1993, A&AS, 97, 355 
Burrows, D. N., Hill, J. E., Nousek, J. A., et al. 2004, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 5165, Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. 
K. A. Flanagan & O. H. W. Siegmund, 201-216 
Camilo, F, Ransom, S. M., Gaensler, B. M., et al. 2006, ApJ, 637, 456 
Carter, J. A. & Read, A. M. 2007, A&A, 464, 1 155 
Davelaar, J., Smith, A., & Becker, R. H. 1986, ApJ, 300, L59 
de Rosa, A., Ubertini, P., Campana, R., et al. 2009, MNRAS, 393, 527 
Gabriel, C, Denby, M., Fyfe, D. J., et al. 2004, in Astronomical Society of the 
Pacific Conference Series, Vol. 314, Astronomical Data Analysis Software 
and Systems (ADASS) XIII, ed. F. Ochsenbein, M. G. Allen, & D. Egret, 
759-+ 

Gabriel, C, Ibarra, A., Gonzalez-Riestra, R., et al. 2008, in XMM-SOC-USR- 
TN, Vol. 17, XMM-SOC-USR-TN, ed. XMM-Newton Science Operation 
Centre and XMM-Newton Science Survey Centre, 1-37 

Garmire, G. P., Bautz, M. W., Ford, P. G., Nousek, J. A., & Ricker, Jr., 
G. R. 2003, in Society of Photo-Optical Instrumentation Engineers (SPIE) 
Conference Series, Vol. 4851, Society of Photo-Optical Instrumentation 
Engineers (SPIE) Conference Series, ed. J. E. Truemper & H. D. Tananbaum, 
28^4 

Gehrels, N., Chincarini, G., Giommi, P., et al. 2004, ApJ, 611, 1005 
Gupta, Y., Mitra, D., Green, D. A., & Acharyya, A. 2005, Current Science, 89, 
853 

Ishida, M., Suzuki, K., & Someya, K. 2007, in JX-ISAS-SUZAKU-MEMO, 

Vol. 11, JX-ISAS-SUZAKU-MEMO, ed. Suzaku Team, 1-7 
Jahoda, K., Swank, J. H, Giles, A. B., et al. 1996, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 2808, Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. 
O. H. Siegmund & M. A. Gummin, 59-70 
Jansen, F, Lumb, D., Altieri, B., et al. 2001, A&A, 365, LI 
Kirsch, M. G., Briel, U. G., Burrows, D., et al. 2005, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 5898, Society 
of Photo-Optical Instrumentation Engineers (SPIE) Conference Series, ed. 
O. H. W. Siegmund, 22-33 
Kokubun, M., Makishima, K., Takahashi, T., et al. 2007, PASJ, 59, 53 
Koyama, K, Tsunemi, H, Dotani, T., et al. 2007, PASJ, 59, 23 
La Palombara, N. & Mereghetti, S. 2002, A&A, 383, 916 
Lebrun, F, Leray, J. P., Lavocat, P., et al. 2003, A&A, 411, L141 
Matheson, H. & Safi-Harb, S. 2010, ArXiv e-prints 
Mitsuda, K, Bautz, M., Inoue, H, et al. 2007, PASJ, 59, 1 
Nevalainen, J., David, L., & Guainazzi, M. 2010, ArXiv e-prints 
Plucinsky, P. P., Haberl, F, Dewey, D., et al. 2008, in Society of Photo-Optical 
Instrumentation Engineers (SPIE) Conference Series, Vol. 7011, Society of 
Photo-Optical Instrumentation Engineers (SPIE) Conference Series 
Slane, P., Chen, Y, Schulz, N. S., et al. 2000, ApJ, 533, L29 
Snowden, S. L. 2002, ArXiv Astrophysics e-prints 
Stephenson, C. B. & Sanduleak, N. 1977, ApJS, 33, 459 
Striider, L., Briel, U., Dennerl, K., et al. 2001, A&A, 365, L18 
Stuhlinger, M., Kirsch, M. G. F, Santos-Lleo, M., et al. 2008, in XMM- 
SOC-CAL-TN, Vol. 52, XMM-SOC-CAL-TN, ed. XMM-Newton Science 
Operation Centre, 1-24 
Takahashi, T., Abe, K., Endo, M., et al. 2007, PASJ, 59, 35 
Tian, W. W. & Leahy, D. A. 2008, MNRAS, 391, L54 
Turner, M. J. L., Abbey, A., Arnaud, M., et al. 2001, A&A, 365, L27 



15 



M. Tsujimoto et al.: X-ray Cross-calibration with G2 1.5-0.9 



Ubertini, P., Lebrun, F., Di Cocco, G., et al. 2003, A&A, 41 1, L131 
Verner, D. A., Ferland, G. J., Korista, K. T., & Yakovlev, D. G. 1996, ApJ, 465, 
487 

Warwick, R. S., Bernard, J., Bocchino, F., et al. 2001, A&A, 365, L248 
Weisskopf, M. C, Brinkman, B., Canizares, C, et al. 2002, PASP, 1 14, 1 
Weisskopf, M. C, Guainazzi, M., Jahoda, K., et al. 2010, ApJ, 713, 912 
Wilms, J., Allen, A., & McCray, R. 2000, ApJ, 542, 914 
Winkler, C, Courvoisier, T., Di Cocco, G., et al. 2003, A&A, 41 1, LI 



16 



